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Numerical simulations of turbulent premixed 
H,/O a /N, flames with complex chemistry 

By M. Baum, 1 T. J. Poinsot, 2 AND D. C. Haworth 3 


Premixed stoichiometric H 2 /0 2 /N 2 flames propagating in twc^dimensional tur- 
bulence have been studied using direct numerical simulation (simulations in which 
all fluid and thermochemical scales axe fully resolved) including realistic chemical 
kinetics and molecular transport. Results axe compared with earlier zer^chemistry 
(flame sheet) and one-step chemistry simulations. Consistent with the simpler mod- 
els the turbulent flame with realistic chemistry aligns preferentially with extensive 
strain rates in the tangent plane and flame curvature probability density functions 
are close to symmetric with near-zero means. By contrast to simple chemistry re- 
sults with non-unity Lewis numbers (ratio of thermal to species diffusivity), local 
flame structure does not correlate with curvature but rather with tangential strain 
rate. Turbulent straining results in substantial thinning of the flame relative to t e 
steady unstrained laminar case. Heat release and H 2 0 2 contours remain thin and 
connected (‘flamelet-like’) while species including H-atom and OH are more dif- 
fuse. Peak OH concentration occurs well behind the peak heat-release zone. This 
work suggests the feasibility of incorporating realistic chemistry into full turbulence 
simulations to address issues such as pollutant formation in hydrocarbon-air flames. 


1. Introduction 

Turbulent premixed combustion in practical devices is a complex phenomenon 
combining chemical kinetics, molecular transport, and hydrodynamic turbulence in 
difficult geometric configurations. Increasing regulative and competitive pressures 
demand improved physical understanding and predictive modeling capability for 
combustion phenomena including ignition, quenching, and pollutant formation, tor 
example, many governments impose limits on emissions of oxides of nitrogen, carbon 
monoxide, and unburned hydrocarbons from gasoline-fueled spark-ignited automo- 
tive engines. Pending regulations in the United States seek further reductions in the 
allowable levels of these pollutants and provide for species differentiation in hydro- 
carbon emissions to account for the relative propensity of different components to 
promote ozone formation in the atmosphere. Many aspects of chemical kinetics (e.g. 
aromatic hydrocarbon chemistry, Westbrook 1991) are poorly understood even at a 
fundamental level and emissions models for engineering applications remain largely 

1 Ecole Centrale, Paris 

2 C.N.R.S., Institut de Mecanique des Fluides de Toulouse, France 

3 General Motors Research & Environmental Staff, Warren, MI 
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empirical. Thus complex chemical kinetics in turbulent premixed combustion is a 
crucial topic for research and modeling. 

Both experimental and numerical investigations of chemical kinetics influences in 
turbulent premixed combustion are complicated by the strong coupling of hydrody- 
namics with thermochemistry and by resolution requirements: hydrodynamic and 
thermochemical spatial and temporal scales span many orders of magnitude in pre- 
mixed flames at high Reynolds and Damkohler numbers (the latter being the ratio 
of a characteristic flow time scale to chemical time scales). Thus direct numerical 
simulation (DNS) of practical turbulent premixed flames is impracticable at present 
and will remain so for the foreseeable future. (Here and in the following, ‘DNS’ is 
taken to mean simulations in which all scales of motion are fully resolved both 
spatially and temporally so that no turbulence modeling — explicit spatial or tem- 
poral filtering — is employed.) Idealized model problems are amenable to numerical 
analysis, and carefully designed numerical simulations are well suited to isolating 
specific phenomena and extracting fundamental physical understanding. However, 
one must remain wary in extrapolating results from idealized model problems to 
practical combustion systems. 

Numerical combustion studies have tended to proceed along one of two parallel 
paths. Either detailed chemical kinetics models have been implemented in simple 
flow configurations (one-dimensional laminar, e.g. Warnatz 1981, Drake & Blint 
1988; axisymmetric laminar, e.g. Smooke et al . 1990), or turbulence simulations 
have been performed with embedded simple-chemistry models (zero-chemistry or 
flame sheet, e.g. Kerstein et al 1988, Girimaji & Pope 1992; zero-heat-release one- 
step Arrhenius, e.g. Rutland et al. 1990; one-step Arrhenius with heat release, e.g. 
Haworth & Poinsot 1992). The detailed chemistry scheme for methane— air oxidation 
employed by Xu & Smooke (1991) includes 26 species and 45 reactions. Solving 
the resulting stiff system of convection/diffusion/reaction equations for a laminar 
axisymmetric Bunsen burner flame requires on the order of 100 CPU hours using 
state-of-the-art numerical algorithms on a high-end workstation (~ 26 Mflops). 
This is the same order of computational resources as that required for DNS of 
nonreacting homogeneous isotropic turbulence at a Taylor-scale Reynolds number 
of Rt\ & 90 (Pope 1991). 

The present research represents an attempt to bridge the two numerical ap- 
proaches by coupling complex chemistry and full turbulence simulations. Here 
^2 — O 2 chemistry (inert N2 diluent) has been modeled using a nine— species, 19- 
reactant scheme (Miller et al 1982) including detailed molecular transport in two- 
dimensional isotropic turbulence. Specific goals are threefold. First, we wish to 
explore the feasibility of complex-chemistry DNS for future application to more 
complex thermochemical systems (e.g. hydrocarbon fuels). Ignition, quenching, 
and pollutant formation issues are of practical interest in such systems. Second, we 
will compare results obtained using the present detailed-chemistry scheme to those 
obtained earlier using simpler zero- and one-step chemistry models. In particular 
we are interested in statistics of strain rate, curvature, and local burning rate along 
the flame front and the robustness of flamelet models with realistic kinetics. Finally, 
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we will compare the local structure of the turbulent premixed flame with detailed 
kinetics to that of an unstrained laminar flame having the same thermochemical 
properties. Questions of interest include possible shifts in radicals relative to lam- 
inar profiles and which radical species might best be suited to marking the region 
of peak heat release in experiments. 


2. Background: DNS of turbulent premixed combustion 

DNS has proved to be a valuable tool in addressing fundamental physical ques- 
tions and in the construction of models for turbulent premixed combustion. The 
past several years have seen DNS applied to a hierarchy of premixed systems of in- 
creasing complexity. Results of these simulations generally have been interpreted in 
the framework of the flamelet regime of combustion wherein a thin laminar-like reac- 
tion zone separates unburnt gas from hot burnt products. (In propagating-interface 
type models, this structure is imposed.) Roughly speaking, flamelet combustion cor- 
responds to conditions where the largest chemical scales are small compared with 
the smallest hydrodynamic scales: however, it has proven to be a worthwhile frame- 
work in which to interpret results even in cases where turbulence microscales are 
smaller than the flame thickness (Haworth & Poinsot 1992). Moreover, flamelet 
models are widely used as combustion ‘submodels’ in Reynolds-averaged computa- 
tions of premixed turbulent combustion in engineering applications (e.g. El Tahry 
1990, Boudier et al. 1992). 

2.1. Propagating-interface (zero-chemistry) models 
Analytic studies of the kinematics of a propagating surface in turbulent flow have 
resulted in relationships that provide a basis for flamelet models and for interpreting 
numerical simulation results. For example, the roles of hydrodynamic straining 
and flame curvature in modifying the area A of a propagating surface element are 
expressed in the relationship (Pope 1988, Candel & Poinsot 1990), 


1 dA Sd 

Ali~ at + n ' 


( 1 ) 


Here at is the hydrodynamic strain rate in the plane tangent to the surface, sj is 
the speed of advance of the propagating surface relative to the fresh gas, and It is 
the radius of curvature (11 < 0 for surface elements propagating towards the center 
of curvature, e.g. concave towards reactants for the premixed flame). This relation- 
ship motivates attempts to isolate strain-rate from curvature effects in turbulent 
premixed combustion, even though the two are not completely independent (Pope 
1988, Haworth k Poinsot 1992). 

Numerical simulations (DNS) incorporating propagating surface models have pro- 
ceeded via a variety of front-tracking algorithms (reviewed in Oran k Boris 1987), 
field-equation approaches (Kerstein et al. 1988, Ashurst et al. 1988), and statistical 
ensembles of infinitesimal surface elements (Girimaji k Pope 1992). These zero- 
chemistry flame-sheet models have contributed significantly to our understanding 
of flame topology and have provided quantitative information on statistical cor- 
relations that is useful in the construction and calibration of models of turbulent 



348 


M. Baum, T. J. Poinsot, & D. C. Haworth 


premixed combustion in the flamelet regime (e.g. Bray & Cant 1991). The influence 
of chemical kinetics effects in such models generally is felt only through the propa- 
gation speed, whose dependence on thermochemical and/or turbulence parameters 
must be prescribed. 

2.2. Simple- chemistry models 

Studies of the dynamic interactions between fluid flow and finite-rate chemistry 
in flamelet or non-fiamelet regimes require that the flame structure be resolved. 
However, the addition of chemical length and time scales that are of the order of or 
smaller than the smallest turbulence scales implies that, for a given spatial and tem- 
poral resolution, a smaller range of hydrodynamic scales can be simulated compared 
to computations in which the flame has no internal structure. Two-dimensional 
vortex methods have been used to study interactions between turbulent fluid flow 
and finite-rate chemistry by a number of authors including Ashurst k Barr (1983), 
Ghoniem k Krishnan (1988), and Ashurst et al. (1987). Three-dimensional simula- 
tions including finite-rate chemistry (constant density, zero heat release, single-step 
Arrhenius chemistry) have been reported by Rutland et al. (1990): results included 
distributions of local burning rate over a range of Damkohler numbers. Further 
three-dimensional constant-density simulations been published by El Tahry et al. 
(1991) Rutland k Trouve (1990) (a study of Lewis number effects) and by Cant et 
al. (1990) (a study of statistics relevant to the Bray— Moss— Libby model of turbulent 
premixed combustion). 

The present work follows a number of two-dimensional simulations with variable 
fluid properties and heat release (Poinsot 1991; Poinsot et al. 1990, 1991, 1992; 
Meneveau k Poinsot 1990, Haworth k Poinsot 1992, Poinsot k Haworth 1992). 
Contributions of these studies include a characterization of the scales of turbulent 
motion that influence flame structure (Poinsot e< al. 1990, 1991), investigations 
of flame quenching (Poinsot et al. 1991, Meneveau k Poinsot 1990), identifica- 
tion of Lewis number effects (Haworth k Poinsot 1992), a study of ignition and 
early flame-kernel growth (Poinsot 1991), and a model for flame-wall interactions 
(Poinsot k Haworth 1992). Compared to three-dimensional constant— property sim- 
ulations, two-dimensional variable— density simulations allow a wider dynamic range 
of scales and full two-way fluid-chemistry coupling. Values of relevant dimension- 
less parameters in the present study are given in Section 3. 

However, the dynamics of two-dimensional turbulence are not identical to those of 
three-dimensional turbulence (Batchelor 1953, Herring et al. 1974, Lesieur 1987). 
In particular, the vortex— stretching mechanism for the cascade of energy to pro- 
gressively smaller scales of motion is absent in two dimensions and the smallest 
scales of motion do not follow the usual Kolmogorov scaling. Thus the statistics of 
small-scale quantities especially are expected to differ between two-dimensional and 
three-dimensional simulations. Partial justification for the appropriateness of two- 
dimensional studies of premixed flame structure can be found in three-dimensional 
results: the topology of a propagating surface in three-dimensional turbulence has 
been found to be primarily two-dimensional, particularly those surface elements 
having the highest curvatures (Ashurst 1990, Cant et al. 1990, Girimaji k Pope 
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1992). That is, a flame tends to be locally cylindrical rather than spherical in 
shape. Further a posteriori justification can be found in the two-dimensional results 
themselves, where many features in common with the three-dimensional findings 
have been noted (Haworth k Poinsot 1992). Common results include: preferential 
alignment of the flame with extensive strain rates in the tangent plane; scaling of 
flame-area-averaged mean tangential strain rate with turbulence micro-timescales; 
and nearly symmetric pdf’s of flame curvature having near-zero mean. 

The present work follows most immediately the study of nonunity Lewis num er 
effects (Le = ratio of thermal to species diffusivity) in two-dimensional turbulent 
premixed combustion reported by Haworth k Poinsot (1992). In addition to the 
results already discussed, it was found that for Le - 1, the local burning velocity of 
the turbulent flame is everywhere nearly identical to that of an undisturbed laminar 
flame; for nonunity Lewis numbers, the local burning velocity differs from that of 
the laminar flame and correlates strongly with the local flame curvature; curvature 
effects cancel out in the mean to leave the mean extensive tangential strain rate as 
the principal influence on the mean burning velocity of the turbulent flame; and 
thermodiffusive effects result in more flame area for Le < 1 than for Le > 1. It was 
further argued that these molecular transport effects should remain important at 
higher Reynolds numbers and for complex chemistry, provided that a single global 
deficient-reactant-based Lewis number can be defined. Here the simulations are 
extended from one-step chemistry with simple transport to realistic chemistry and 
transport, permitting a direct assessment of the second claim. We continue to adopt 
a flamelet viewpoint for diagnostics and analysis. 

3. Problem definition 

8.1. Governing equations and numerical methods 

The set of equations solved is the compressible multi-species reacting flow equa- 
tions comprising conservation of mass, linear momentum, energy, and N s species 
mass fractions. In Cartesian tensor notation (no sum on Greek indices), 


dp dpuj _ 

dt + dxj 

( 2 ) 

dpui dpUiUj dp drij 

dt + dxj dxi dxj 

( 3 ) 

dpe t d(pe t + p)uj d(uj T ij) fob 

dt 4 dx } dxi dxj ' 

( 4 ) 

dpY a dpY 0 n } dpY a V a j w 

dt dxj dxj 

( 5 ) 


Here u , is the i th component of the fluid velocity, et is the total energy density per 
unit mass, and Y a , a = 1, . . . , N„ is the mass fraction of species a. The fluid mass 
density is p, p is the thermodynamic pressure, and r tJ and q } are, respectively, the 
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viscous stress tensor and the heat flux vector. The molar chemical production rate 
of species a is u> a while W a is its molecular weight. The variable V a j is the diffusion 
velocity for species a. 

The set of governing equations is closed with the ideal gas equation of state and 
constitutive relations, 


P = pRT , 

(6) 


_ f dui duj 2 . du*\ 



T>} ~ 11 [dx, + dxi ~ 3 6ij dT k ) ' 

(7) 


1 



ei = 2 U * U * + e i 

(8) 

- 

N, 


- 

e = T h a Y a - ? , 
a=l P 

(9) 

1 

f Tl 


r 

h a (T 1 ) = Ah° fa + C pa (T)dT , 

Jt 0 

(10) 

m 

_ 

Q) ^ Qx ^ ^ I haYaVaj . 

3 a=l 

(11) 



In these equations, JI and A are, respectively, the mean viscosity and thermal conduc- 
tivity. The formation enthalpy of species a at the reference temperature To is Ah°, a , 
and h a is the enthalpy of species a. (Here the reference state is To =298.15 K.) 

Species molecular transport is modeled using Fick’s law. An effective species 
diffusivity D a is computed as follows: 


D a 


( 1 zM 

N, 

X) Xp/D a0 


( 12 ) 


where D af} are the binary diffusion coefficients and X a is the mole fraction of species 
a. Diffusion coefficients are related to the diffusion velocity V aj by, 


— 



(13) 


The diffusion velocities resulting from Eqs. (12) and (13) do not, in general, satisfy 
conservation of mass. A correction velocity V cj is added such that the net diffusive 
flux is equal to zero, 

N. 

Y,Y a V a j = 0 . 

a=l 


( 14 ) 
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periodic 




FIGURE 1. Schematic of computational domain: a premixed flame separates 
unburnt reactants from burnt products. 

Finally, we have, 


V Qj = V Q , + V cj = 


'aj 


c] 


E Yp(Dp/Xp)Xpj 
D a Y 

- ^ N. 

E Yf 

0=1 


(15) 


Species production rates are given by the Arrhenius law with (forward) rate 
constants kf a of the form, 

kfa = A a T exp(-E 0a /RoT) , (16) 

where Ro is the universal gas constant. Hydrogen-oxygen kinetics has been mod- 
eled using the nine-species, 19-reaction scheme devised by Miller et al (1982). 
Coefficients for this scheme are summarized in Table I. Fluid properties, molecular 
transport coefficients, and reaction source terms are computed using CHEMKIN 
and TRANSPORT (Kee et al. 1980, Kee et al. 1983). 

Using these assumptions and a Cartesian frame of reference, the conservation 
equations are solved using a high-order finite-difference scheme (Lele 1992). The 
calculations are initialized with reactants on one side of the computational domain 
and products on the other; these are separated by a laminar premixed flame (Figure 
1). The initial laminar profiles are themselves steady one-dimensional solutions to 
Eqs. (2)-(16). On lateral boundaries, periodic conditions are enforced while non- 
reflecting boundary conditions are used on inflow/outflow boundaries (Poinsot & 
Lele 1992). Isotropic two-dimensional turbulence is prescribed with a turbulence 
spectrum that is the same as that used in earlier studies (e.g. Haworth &; Poinsot 
1992). Two parameters suffice to define the initial energy spectrum: the rms turbu- 
lence velocity u p = u'(t = 0) and the peak energy wavelength L,. These and other 
relevant parameters are summarized in Table II. 
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Table I. Reaction mechanism rate coefficients in the form k f = AT P exp{- E 0 / RoT) 
(Eq. 16) (Miller et al. 1982). Units are moles, cubic centimeters, seconds, degrees 
K, and calories/mole. Species: H 2 0 2 OH O H H 2 0 H0 2 H 2 0 2 N 2 . 


# 

Reaction 

A 

0 

E 0 

Ri 

H 2 + 0 2 ^ 2 OH 

1.7 x 10 13 

0 . 

47780 

r 2 

H 2 + OH ^ H 2 0 + H 

1.17 x 10 9 

1.3 

3626 

r 3 

H + 0 2 ^ OH + 0 

5.13 x 10 16 

-0.816 

16507 

R 4 

0 + H 2 ^ OH + H 

1.8 x 10 *° 

1.0 

8826 

r 5 

H + 0 2 + M ^ H0 2 + M“ 

2.1 x 10 1 ® 

- 1.0 

0 

Re 

H + 2 O 2 ▼=* HO 2 + O 2 

6.7 x 10 19 

-1.42 

0 

Rt 

H + 0 2 + N 2 ^ H0 2 + N 2 

6.7 x 10 19 

-1.42 

0 

Rs 

OH + H0 2 ^ H 2 0 +0 2 

5 x 10 13 

0 . 

1000 

R 9 

H + H0 2 ^ 2 OH 

2.5 x 10 14 

0 . 

1900 

Rio 

0 + H0 2 ^ 0 2 + OH 

4.8 x 10 13 

0 . 

1000 

Rn 

2 OH ^ 0 + H 2 0 

6 x 10 ® 

1.3 

0 

Rl 2 

h 2 + m^h + h-i-m 6 

2.23 x 10 12 

0.5 

92600 

Rl3 

o 2 + m^o + o + m 

1.85 x 10" 

0.5 

95560 

Rl4 

H + OH + M ^ H 2 0 + M c 

7.5 x 10 23 

- 2.6 

0 

Rl5 

H0 2 + H H 2 + 0 2 

2.5 x 10 13 

0 . 

700 

Rl6 

2 H0 2 H 2 0 2 + 0 2 

2 x 10 12 

0 . 

0 

R]7 

H 2 0 2 + M ^ OH + OH + M 

1.3 x 10 17 

0 . 

45500 

Rl8 

H 2 0 2 + H ^ h 2 + ho 2 

1.6 x 10 12 

0 . 

3800 

Rl9 

h 2 o 2 + oh ^ h 2 o + ho 2 

1.0 x 10 13 

0 . 

1800 


“ Third-body efficiencies: fc(H 2 ) = 3.3fc(Ar), ifc(H 2 0) = 21.0Jfc(Ar). 

6 Third-body efficiencies: £(H 2 ) = 3.0fc(Ar), fc(H) = 2.0ifc(Ar), ifc(H 2 0) = 6.0fc(Ar). 
c Third-body efficiencies: &(H 2 0) = 20.0Jb(Ar). 

In Table II, a subscript V refers to properties in the unburnt reactants. Key 
parameters are: the rms turbulence intensity u ' , the turbulence integral length 
scale / based on two-point velocity correlations and timescale r = l/u', and flame or 
chemical scales (the unstrained laminar flame speed), S, = (T b -T u )/(dT/dx\ mai ) 
(the laminar flame thickness), and tj = Si / s® . The fuel mass equivalence ratio is $. 
There are n x x n v computational grid points in the L boz x L box square computational 
domain. 

The turbulence Reynolds numbers and ratio of rms turbulence intensity to lami- 
nar flame speed of Table II imply a regime of combustion where turbulence is intense 
compared to the laminar flame speed, and flames are thinner than the turbulence 
integral scale but are thicker than turbulence microscales. Turbulence and flame 
timescales are of the same order. 
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Table II. Initial parameters for two-dimensional simulations. 


Case 

$ 

T u 

P 

-? 

6i 

Re Li 

Rei 

«'/*? 

LJ6, 

1/6, 

r/r; 

Lbox 

n x 






tip L, 







71 y 


H 

[K] 

[kPa] 

[m/s] 

[cm] 

[-1 

H 

H 

H 

H 

(-] 

[cm] 

[-] 

i 

1.0 

800 

101 

19.5 

0.044 

2000 

650 

2.1 

6.3 

2.0 

0.48 

2.00 

487 

2 

1.0 

800 

101 

19.5 

0.044 

4000 

1300 

4.2 

12.7 

4.0 

0.95 

1.24 

301 


S.2. Diagnostics 

Postprocessing of the two-dimensional computed fields (snapshots at fixed times) 
begins by defining a flame front as an isocontour of reactant mass fraction or tem- 
perature. Here we have used Y 0 , = 0.22, corresponding to the position of peak 
heat release in the laminar flame (Figure 2). Once the flame front has been located, 
the local normal and local curvature are readily computed. Curvatures concave 
towards the hot products are taken to be positive. One-dimensional cuts normal 
to the flame are taken: it is these profiles that define the local ‘structure’ of the 
turbulent flame. We compare the local turbulent flame profiles with the steady 
unstrained one-dimensional laminar flame profile for the same chemistry and fluid 
properties. Of particular interest is the distribution along the flame of the normal- 
ized local flame speed (‘flamelet speed’) s„ defined as the integral of the heat-release 
profile in a direction locally normal to the flame. The local heat-release rate wq is 
defined by, 

Wq ~ WaAhfa > (^) 

0=1 

and the normalized flamelet speed is then, 

s n = J dn / s„o • (1®) 

A normalized flame length C* and turbulent flame speed s* T are defined as, 

C = L flame/ L flame 0 i (^) 

Sj- = WQ / WQO = (•*!»)■£* • (20) 

In Eqs. (18)-(20), the subscript ‘0’ denotes the value for the unstrained planar 
laminar flame; overbars denote volume (area) averages while angled brackets are 
reserved for flame-area-weighted means (length-weighted in two spatial dimen- 
sions). A check of self-consistency is to verify the equality between s * T calculated 
in the two ways given in Eq. (20): they are found to agree to within a few percent 
in all cases. 
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Figure 2. Normalized (by their maxiumum value) steady unstrained laminar 
flame profiles: a)-c) complex chemistry (Tables I and II); d) one-step chemistry 

w = ApVexp( 0 = (T - T u )/(T h - T„), a = 0.75, P = 8.0, A = 146 

(Haworth & Poinsot 1992). 

4» Results 


4-1- Laminar profiles 

Profiles of the steady one-dimensional laminar solution are displayed in Fig- 
ure 2a-c, where each quantity has been normalized by its maximum value. Cor- 
responding profiles for the one-step Arrhenius scheme used in Haworth & Poinsot 
(1992) are shown in Figure 2d. The abscissa is in units of the laminar flame thick- 
ness 8 /, which is of the order of the heat-release zone thickness. A striking difference 
between the simple- and complex-chemistry temperature and species profiles is the 
extended tails into the burnt gas for the complex chemistry case. One might an- 
ticipate difficulty in maintaining this structure in turbulent flames: even in cases 
where all turbulence scales are larger than 8 { , turbulence may still be able to modify 
the laminar flame structure on the burnt side. Other noteworthy features of these 
profiles include the significant ; leakage’ of fuel and oxidizer into the products, high 







H 2 / 0 2 /N 2 flame. 


355 



FIGURE 3. Computed fields at one instant of time, Case 1, time t/rj — 1.24: 
a) vorticity contours with flame isocontour in magenta; h) temperature; c) heat 
release; d) H2O2 mass fraction; e) H mass fraction; f) OH mass fraction. Contour 
color scale ranges from blue (min) to red (max). In a), the full computational box 
is shown while b)-f) depict only the subdomain indicated by the black outline in 
a). 
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OH concentrations in the burnt gas, and locations of peak radical concentrations 
relative to peak heat release. The high product temperature (2840 K) is responsible 
for the high equilibrium H 2 , 0 2 , and OH concentrations on the burnt side. Peak OH 
concentration occurs well behind the peak heat— release zone, while H— and O atom, 
H 2 0 2 , and H0 2 peak closer to the location of peak heat release. 

4.2. Turbulent flame contours 

Examples of computed fields at one instant of time for Case 1 at t/rf = 1.24 
are given in Figure 3 (unburnt gases are on the left, burnt products are on the 
right). The initially planar flame has been severely distorted by the turbulence by 
this time, even to the extent of tearing off pockets of fresh gas which are carried 
into the product side. Eventually these pockets of fresh gas are consumed for 
these adiabatic flames. Nevertheless, the structure of the underlying laminar flame 
remains readily recognizable. Heat release contours (Fig 3c) suggest a thin flamelet— 
like structure albeit with some islands of fresh gas burning enveloped by the hot 
products. Some radicals including H 2 0 2 (Figure 3d) show similar behavior to the 
heat release: profiles remain thin and appear to mark the region of maximum heat 
release. On the other hand, H-atom levels (Figure 3e), while peaking close to the 
maximum heat-release zone, show a long diffuse tail on the burnt side, presumably 
a consequence of the high H diffusivity. Gradients in OH mass fraction are steep 
through the region of maximum heat release (Figure 3f) but peak OH occurs well 
behind the flame front and isocontours of high OH concentration do not remain 
connected on the burnt side. Most of these features might be anticipated from 
Figure 2a-c. Thus in the present case, OH does not serve as a good marker of the 
region of maximum heat release in the flame. 

An example of the influence of turbulence parameters on flame structure is given 
in Figure 4. There the isocontour used to define the flame at one instant of time is 
plotted for Cases 1 and 2 at equal flame times t/rf = 0.36 and at equal turbulence 
times t/r 0 « 0.8 (r 0 = r(t = 0)). The smaller integral scale of Case 1 yields a. more 
highly wrinkled flame that, as illustrated in Figure 3, even shows pockets in the 
burnt gas at later times. Case 2, by contrast, yields a flame with large-scale folding 
versus fine-scale wrinkling. 

4.3. Statistics of flame curvature, strain rate, and flamelet speed 

Statistics of tangential strain rate, flame curvature, and local flamelet speed are 
presented in Figs. 5-9. Probability density functions (pdf’s) of flame tangential 
strain rate and flame curvature are given in Figs. 5 and 6, respectively. Figure 5 
confirms the findings of several earlier studies using zero- and single-step chemistry 
models: the turbulent flame aligns preferentially with extensive strain rates in the 
tangent plane (Ashurst 1990, Girimaji k Pope 1992, El Tahry et al. 1991, Cant 
et al. 1990, Rutland et al. 1990). Figure 6 also is reminiscent of results found 
using simpler models (Haworth & Poinsot 1992) and experiments (Lee et al. 1991): 
curvature pdf’s are close to symmetric with near-zero mean curvature. While some 
skewing of the pdf’s towards negative curvatures might be expected at low u'/s'j 
(Rutland & Trouve 1990, Becker et al. 1990), no such skewing is evident here. The 
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a) t/r F = 0.36 b) t/r 0 = 0.8 




Figure 4. Isocontours Yo 2 = 0.22 for Cases 1 and 2 in a 1.24 cm x 1.24 cm 
region: a) t/rf = 0.36; b) t/ro « 0.8. 

more highly wrinkled nature of the Case 1 flame is manifested in the broader tails 
(higher rms) of the curvature pdf for Case 1 versus Case 2. 

Scatter plots of local flamelet speed s n versus local tangential strain rate and 
local flame curvature are given in Figs. 7 and 8, respectively. Each point in these 
two figures corresponds to a point sampled uniformly along the isocontour defining 
the turbulent flame. While no correlation of local burning velocity with curvature 
is evident, the decrease in s n with increasing tangential strain rate is surprisingly 
strong. In fact, this correlation develops over time in the computations: at time 
t = 0, the normalized flamelet speed s n is everywhere equal to unity. The flame- 
area-averaged normalized mean flamelet speed (s„) decreases with time (Figure 9). 
Figure 9b, for example, portrays the evolution from the initial delta function pdf of 
flamelet speed at time t = 0 for Case 2: both the mean value ($„) and the value of 
s n at which the pdf peaks decrease with time as the flame develops (Table III). 

No statistically steady state is achieved in the duration of the present compu- 
tations. While the mean consumption rate per unit area of flame (s n ) (Eq. 18) 
is significantly smaller than unity, the increase in flame area C * (Eq. 19) more 
than compensates to yield normalized turbulent flame speeds sf (Eq. 20) that are 
greater than unity, as expected. The evolution of these global quantities with time 
is summarized in Table III. 

Peak heat-release rates (not shown) remain within ±5% of the laminar value in 
all cases, suggesting that this reduction of local consumption rate ( s n ) is primarily 
a flame thickness effect. These results are unexpected in view of the findings in Ha- 
worth & Poinsot (1992) for simple chemistry. There it was found that, for nonunity 
Lewis numbers, curvature dominated the local flame structure, and a correlation 
between local burning velocity and local tangential strain rate was evident only for 
Le — 1. In no case was the straining sufficiently high to yield the extremely low val- 
ues of (s„) found here. While the fuel-based (H 2 ) Lewis number in the present case 
is clearly less than unity, no single global deficient-reactant-based Lewis number 
can be unambiguously defined for $ = 1 (Chelliah & Williams 1987). Arguments by 
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tang, strain rate x flame time 

FIGURE 5. Pdf of normalized tangential strain rate, Case 1, t/r / = 0.36: 
mean=1.816, rms=1.394. 



FIGURE 6. Pdf’s of normalized flame curvature, t/ro « 0.8. Case 1: mean=0.0056, 
rms=0.82; Case 2: mean=0.0043, rms=0.35. 


Joulin & Mitani (1981) suggest an equivalent Lewis number of somewhat less than 
unity for the present thermochemistry: however, the decreasing flamelet speed with 
increasing tangential strain rate of Figure 7 suggests that, if anything, the effective 
Lewis number may be greater than or equal to unity here. 




360 


M. Baum , T. J. Poinsot , & D . C. Haworth 



fang, strain rah x flams firm fang, strain rah x flams fkns 


FIGURE 7. Scatter plots of normalized flamelet speed s n (Eq. 18) versus normalized 
tangential strain rate at t/ro ~ 0.8: a) Case 1; b) Case 2. 


Table III. Evolution with time of normalized flamelet speed { s n ) (Eq. 18), normal- 
ized flame area C* (Eq. 18), and normalized turbulent flame speed s ^ (Eq. 20) for 
Case 2. 


; t/r 0 

(Sn) 

C* 

Sp 

i 

0.0 

1.00 

1.00 

1.00 

0.23 

0.96 

1.10 

1.06 

j 0.36 

0.91 

1.24 

1.13 

| 0.80 

0.75 

1.86 

1.40 


4-4- Local flame structure 

Profiles of heat-release rate and OH mass fraction along the local flame normal 
are overlaid on the corresponding laminar profiles in Figure 10 for Case 2 at two 
instants of time. The turbulent profiles have been uniformly sampled along the 
turbulent flame front; only nine or ten turbulent profiles are shown, for clarity. 
At the earlier time (t/r 0 = 0.23), local turbulent profiles collapse neatly onto the 
steady unstrained laminar contours. At later times ( t/r 0 = 0.80), the collapse 
remains reasonable on the fresh-gas side while the turbulent profiles generally are 
shifted inwards towards the peak heat-release zone on the burnt-gas side. This 
appears to be a consequence of the net extensive tangential strain rates thinning 
the flame in the turbulent flow and reducing the local flamelet speed while peak 
heat release rate remains relatively unchanged. 
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FIGURE 8. Scatter plots of normalized flamelet speed s n (Eq. 18) versus normalized 
flame curvature at </tq « 0.8: a) Case 1; b) Case 2. 



FIGURE 9. Pdf’s of normalized flamelet speed s n (Eq. 18): a) Case 1; b) Case 2. 


5. Discussion 

Predictive models of turbulent premixed combustion for engineering applications 
remain at an early stage of development. Key physical processes including ignition 
and early flame kernel development, complex hydrocarbon- air reaction kinetics, 
turbulence-chemistry interactions, and flame-wall interaction remain poorly un- 
derstood even at a fundamental level. Direct numerical simulation has played an 
increasingly important role in recent years both in contributing to fundamental un- 
derstanding and in providing guidance for modeling. Fundamental contributions 
of DNS include elucidation of flame topology and structure (Kerstein et al 1988, 
Ashurst et al. 1988, Rutland k Trouve 1990, Ashurst 1990, Cant et al. 1990, 
Rutland et al. 1990, El Tahry et al. 1991, Girimaji k Pope 1992), a characteri- 
zation of the scales of turbulent motion that influence flame structure (Poinsot et 
al 1990, 1991), flame quenching (Poinsot et al 1991, Meneveau k Poinsot 1990), 
and thermodiffusive (Lewis number) effects (Rutland k Trouve 1990, Haworth k 
Poinsot 1992). Moreover, DNS results have been used to construct and to cali- 
brate models for premixed flame propagation (Bray k Cant 1992), ignition and 
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FIGURE 10. Unnormalized laminar and turbulent flame profiles for Case 2: a) 
heat-release rate, f/r 0 = 0.23; b) heat-release rate, f/r 0 = 0.80; c) OH mass fraction, 
</tq= 0.23; d) OH mass fraction, t/r 0 = 0.80. 


early flame-kernel growth (Poinsot 1991), and flame-wall interactions (Poinsot & 
Haworth 1992). These models can be incorporated into multidimensional Reynolds- 
averaged modeling approaches for engineering devices (e.g. Boudier et al 1992). 

The present simulations represent a first attempt to incorporate realistic chemical 
kinetics and transport into full turbulence simulations. Principle results are: 1) 
Realistic chemical kinetics and transport have been implemented successfully into 
full turbulence simulations. 2) Consistent with the simpler models, the complex- 
chemistry turbulent flame aligns preferentially with extensive strain rates in the 
tangent plane and flame curvature pdf’s are close to symmetric with near-zero 
means. 3) By contrast to simple chemistry results with non-unity Lewis numbers, 
local flame structure does not correlate with local flame curvature but rather with 
local tangential strain rate. Turbulent straining results in substantial thinning of the 
flame relative to the steady unstrained laminar case. Peak heat-release rate remains 
close to the laminar value, but integrated heat release along the local flame normal 
is substantially lower in the turbulent flame. 4) Heat release and H 2 0 2 contours 
remain thin and connected (‘flamelet-Iike’) while species including H-atom and OH 
are more diffuse. In particular, peak OH concentration occurs well behind the peak 
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heat-release zone. 

The strong correlation of s n with strain rate rather than with curvature (Figures 
7 and 8) for the present computations suggests an effective Lewis number that is 
greater than or equal to unity. However, no amount of straining of the simple- 
chemistry flames in Haworth k Poinsot (1992) resulted in the extreme thinning and 
decrease in (s n ) than has been observed here. For the present stoichiometric reac- 
tants, definition of a single deficient-reactant-based Lewis number is problematic 
(Joulin k Mitani 1981, Chelliah k Williams 1987). Computations with fuel-lean 
reactants may help to shed light on these results. 

Present turbulence Reynolds numbers and the normalized turbulence intensity 
u'/s? (Table II) are typical of values found, for example, in reciprocating internal 
combustion (IC) engines at low engine speeds. The computed flames remain thick 
compared to the l j 6 \ ~ 30 typical of propane— air mixtures under standard IC— 
engine operating conditions (stoichiometric, p = 5 atm, T u = 600 K, undiluted). 
However, I /61 can be as low as three for high T u , low p, or high exhaust-gas dilution 
engine cases (Mantzaras et ai 1988; Blint, 1988, 1990). The present choice of 
initial turbulence spectrum (Haworth k Poinsot 1992) results in too little energy 
at low wavenumbers compared with classic equilibrium turbulence: thus there is 
little dissipation in the hot post-flame gas in spite of the high temperature-induced 
viscosity there. Computations with a more realistic initial spectrum are planned. 

The high adiabatic flame temperature (2840 K) for the present thermochemistry 
results in significant equilibrium levels of fuel and oxidizer in the burnt gas and in 
high equilibrium radical concentrations of OH, H-atom, and O-atom (Figure 2). 
Thus it is not surprising that OH serves as a poor marker of the heat-release zone 
for this flame. The H 2 0 2 radical serves better in this regard computationally, but 
may be a poor choice in practice: peak H 2 0 2 mass fraction in the present case 
is 0.000094 versus 0.0348 for OH. Moreover, the reaction-rate constants involving 
H 2 0 2 in Table I are not well established. Nonstoichiometric and cooler reactants 
would result in lower temperatures and narrower OH radical profiles. In these 
cases, OH may better serve as a marker of the reaction zone consistent with the 
experimental observations of Becker et al (1990), for example. 

Further hydrogen-air work thus will include modifications to the initial turbulence 
spectrum, fuel— lean cases, runs to later times t/rj, and additional postprocessing to 
look at relative shifts in species mass fraction profiles relative to the steady laminar 
flame. Additional chemistry for NO x production and extensions to pollutant forma- 
tion in hydrocarbon-air systems are anticipated subjects for future investigations. 
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